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Abstract 

In this paper, we present the problem of counting magic squares 
and we focus on the case of multiphcative magic squares of order 4. 
We give the exact number of normal multiplicative magic squares of 
order 4 with an original and complete proof, pointing out the role 
of the action of the symmetric group. Moreover, we provide a new 
representation for magic squares of order 4. Such representation allows 
the construction of magic squares in a very simple way, using essentially 
only five particular 4x4 matrices . 



1 Introduction 

A magic square is defined as an n x n matrix of integers, where the sum 
of the numbers in each line (i.e., in each row, in each column and in each 
diagonal) is the same. Magic squares have a very rich history (see, e.g., 
the beautiful book of Descombes [Bj) and they have many generalizations. 
There are a lot of contributes to their theory from several different fields 
(see, e.g., the classical book of Andrews [T]). 

A very difficult job is counting how many magic squares there are for a given 
order and a given line sum, but if we relax the requirements and we ask for 
the so-called semi-magic squares (i.e., magic squares without the condition 
on the diagonals), the counting is easier. For order 3, Bona [2] has found 
that the number of semi-magic squares is 

r + 4\ /r + 3\ /r + 2 

4 )n 4 )n 4 

where r is the line sum (let us observe that in [2] semi-magic squares are 
called magic squares). For a general theory of counting semi-magic squares 
you can see [3j, chapter 9. 
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The classic magic squares, which are called normal, are nxn matrices whose 
entries consist of the numbers 0, 1, n-^ — 1 and in each line the sum of the 
numbers is the same. Their number is known for the orders 3, 4 and 5. The 
case of order 4 is especially interesting. There are exactly 880 normal magic 
squares of order 4 (up to symmetries of the square, i.e., in total they are 
880 X 8 = 7040). They were enumerated for the first time in 1693 by Frenicle 
de Bessy. The Frenicle method was analytically expanded and completed 
by Bondi and Ollerenshaw [4]. Of course, one can think to different kinds 
of magic squares, only changing the operation. Then a multiplicative magic 
square is an n x n matrix of integers, where the product of the numbers in 
each line is the same. This is not a new idea, for example multiplicative 
semi-magic squares of order 3 are studied in |7|. There are several differ- 
ent methods to construct multiplicative magic squares. The most obvious 
method is to transform an additive magic square A = (aij) into M = (b""'^), 
for any base b. More intriguing ways are explained in [6] (multiplicative 
magic squares are usually constructed by using orthogonal latin squares or 
geometric progressions) . 

In [8|, the problem of finding the exact number of multiplicative magic 
squares of order 4 has been posed. In order to solve such question, in [9] 
a correspondence between additive and multiplicative magic squares is pro- 
posed. The authors claimed that such function is an isomorphism even 
though it is only an injective morphism, as pointed out in [TO]. Since such 
revision is only a note, the above problem is not completely clarified and the 
correct number of multiplicative magic squares is given without a rigorous 
proof. In the next section, we clarify the above questions, proposing an orig- 
inal and complete proof. Moreover, in the last section a novel representation 
for some magic squares is proposed, allowing an easy way to construct them. 



2 Counting multiplicative magic squares of order 
4 

A normal multiplicative magic square nxn is created using n primes pi, ■■■,Pn- 
All the divisors k = pi - ■ ■ pn must appear once time in the matrix and in 
each line the product is the same. The magic constant is k'^ and it is well- 
known that in a normal multiplicative magic square, in each line, any prime 
compares with power one and exactly two times. For example, a normal 
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multiplicative magic square of order 4 is 

/ 1 P3 • P4 PI-P2 P1-P2-P3- P4\ 

P1-P3- PA P1-P2- P3 PA P2 

P1-P2- Pa Pi P2 ■ P3 ■ Pa P3 
\ P2-P3 P2- Pa Pi ■ P3 Pi ■ Pa ) 

In the following, we use 'M.k,n to indicate the set of normal multiplicative 
magic squares n x n with magic constant k"^. An hard problem to solve 
is to find the order of Mfc^4, as posed in |8j. In order to answer to this 
question, M^^a is made in correspondence with the set of normal additive 
magic squares 4 x 4 [S]. 

A normal additive magic square is a matrix n x n, whose entries consist of 
the numbers 0, 1, ...,n^ — 1 and in each line the sum of the numbers is the 
same. The magic constant is "^""^ ■ Prom now on, we name A„ the set of 
normal additive magic squares n x n. 

It is well-known that the order of A4 is 7040, or 880 up to symmetries of 
the square [3]. Can we find a similar result for Mfc 4? 
In [9] the following correspondence is defined: 

/ : Mfc,„ A2„/2. 

If / is bijective, |Mfc^4| is determined. The function / maps M = (rriij) € 
M.k^n into A = (oij) G A2n/2. The element ruij is made in correspondence 
with Oij, where fljj is the base 10 number , between and — 1, whose base 
2 expression is the n-string associated to ruij. This string is constructed as 
follows: arrange the n prime numbers in ascending order; if pk is a factor of 
ruij place a 1 in the A:th position of the string ; if not, place a 0. 
Such a function / is surely injective, but unfortunately it is not surjective 
and the order of 4 is not IA4I = 7040. Indeed, we can consider the group 
S4 which permutes the four primes involved in WI^^a- 

Theorem 1. The order of divides the order ofW.k,A- 

Proof. Surely when a permutation p £ S4 acts on M G Mfc^4, then p{M) is 
still an element in Mfc^4. Furthermore, if p is not the identity, it does not fix 
any element of 4, i.e., any orbit has 24 elements. Therefore by Burnside 
Lemma |S'4| must be a divisor of |Mfe^4|. □ 

As immediate consequence of the previous Theorem, 7040 can not be the 
order of Mfc^4, since 24 does not divide 7040. The correct answer is given in 
the following theorem. 
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Theorem 2. The order o/Mfc,4 is 4224. 



Proof. If we consider ^ G A4 and we write its entries in base 4, then in each 
hne the sum of the digits in the first position multiphed by 4 and added to 
the sum of the digits in the second position must yield 30. The only possible 
combinations are 

30 = 7x4 + 2 = 6x4 + 6 = 5x4 + 10. 

So if A has a line in which the digits in the second position have sum 2, we 
are in the only two situations 

•0 -0 -1 -1 
•0 -0 -0 -2 

where • indicates the first position of the number in base 4. Considering 
the representation in base 2 of these strings, we have one of the following 
situations 

• • 00 • -00 • -01 • -01 
••00 ••OO ••OO ••lO . 

Now if we try to create M G M;i.,4) using the inverse of /, if A has a line which 
presents one of these situations we can not have a normal multiplicative 
magic square. In fact, if we are in the first situation, for example, then the 
prime P3 will not appear in this line. Similarly if a line of A has digits in 
the second position with sum 10, the possible cases are 

■3 -3 -3 1 
■2 -2 -3 3 

which correspond in base 2 to the strings 

••11 ••ll ••ll ••Ol 
••10 ••lO ••ll ••ll . 
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Once more, if we try, we fail to obtain M G Mfc^4 . 

Finally, if A has a line whose digits in the second position have sum 6, we 
have the possibilities 



•0 


• 1 


• 2 


CO 


•0 


•0 


• 3 


• 3 


•1 


• 1 


• 2 


• 2 


•0 


• 2 


• 2 


• 2 


•1 


• 1 


• 1 


• 3 



With similar arguments as the ones used before, it is easy to see that the last 
two situations do not allow to obtain a normal multiplicative magic square. 
Therefore the squares ^ G A4 from which we obtain a normal multiplicative 
magic square are only those with lines whose entries in base 4 have digits in 
the first and second position composed by strings 



•0 


• 1 


• 2 


• 3 


•0 


•0 


CO 


• 3 


•1 


• 1 


• 2 


• 2 



In [3], all the squares in A4 are classified. The squares in A4 generated 
using only these strings, i.e., which correspond to squares in Mfc_4, are only 
those in category one ([3], p. 510). They are exactly 528 unless of symmetries 
of the square. Therefore our arguments and the injectivity of / allow us to 
conclude that the order of Mfc_4 is 528 x 8 = 4224. □ 

3 A new representation for magic squares 

In this section we see a representation for normal additive magic squares 
which correspond to normal multiplicative ones. 

In [5] a lower bound for the distance between the maximal and minimal 
element in a multiplicative magic square is given. In order to do that, the 
relation between additive and multiplicative magic squares is highlighted, 
recalling that a multiplicative magic square can be factorized as IliPf') 
where Aj's are additive magic squares. Moreover, focusing on magic squares 
of order 4, in [5] the Hilbert basis (composed by 20 magic squares) for such 
magic squares is explicited. However, such representation and the relation 
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between additive and multiplicative magic squares are not really manageable 
in order to construct additive and multiplicative magic squares of order 4. 
Here, the proposed representation allows to determine all (and not only) 
the normal additive magic squares of order 4, which corresponds to normal 
multiplicative magic squares. In this way, they can all be easily constructed 
essentially using only 5 basic matrices. 

Wc consider M = {rriij) € Mfc^4, by means of / we have the correspondent 
A = {aij) G A4 and we consider its entries in base 2. Now we decompose A 
into four matrices ^1,^2,^3,^4, whose entries are only or 1, so that the 
entries in position ij of the matrices ^1, A2, ^3, A4 form the string Oy. 

Example 1. Prom the normal multiplicative magic square 

/PlP2P3 PSPA 1 PlP2P4\ 

Pi P1P3P4: P2P3 

PlPi P2 P2P3Pi PlP3 ' 

V P3 P1P2P3PA PIP2 PA J 

using /, we obtain the normal additive magic square 



/lllO 


0011 


0000 


1101\ 




/14 


3 





13\ 


0101 


1000 


1011 


0110 




5 


8 


11 


6 


1001 


0100 


0111 


1010 




9 


4 


7 


10 


\0010 


1111 


1100 


0001/ 




V2 


15 


12 


1/ 



and it can be decomposed into 



8 



A 








i\ 




A 








A 




n 


1 











1 





i\ 





1 


1 





+ 4 


1 








1 


+ 2 








1 


1 


+ 


1 





1 





1 








1 





1 


1 











1 


1 


1 





1 





VO 


1 


1 


0) 




^0 


1 


1 


0) 






1 





0) 




^0 


1 





1/ 



Since A is derived from a matrix in Mjt^4, by means of /, the matrices 
Ai,A2,As,A4 are all and only those having in each line exactly two ones. 
We call forms these matrices which allow us to construct any magic square 
in A4 corresponding to magic squares in Mj^ 4. These forms are exactly 16. 

Theorem 3. There are 16 different matrices 4x4, with entries or 1, such 
that in each line there are exactly two ones. 

Proof. The strings that we can use to generate these forms are only six: 
1100, 1010, 1001, 0110, 0101, 0011 . 
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If we choose as a diagonal a string with different extremes, then we have 
only two possible different forms. For example: 

/I • • A /I • • 0\ /I 1 0\ 
1 • • • 1 • • 10 1 

•0-^--0-^0101 

\1 1 



/I 



0/ 

A 

0/ 



n 



0/ 

0/ 



/I 



1 



0/ 

1\ 



1 

0/ 



On the other hand, if we choose as a diagonal a string with same extremes, 
then we have four possible forms. Since we have four strings with different 
extremes and two strings with same extremes, our forms are 4 • 2 + 2 • 4 = 
16 . □ 

Prom these 16 forms we can individuate five fundamental forms: 



/o 





1 


1\ 





1 





1 


1 





1 





VI 


1 





0/ 


/o 





1 


1\ 


1 


1 








1 


1 








VO 





1 


1/ 


/o 


1 





1\ 


1 


1 














1 


1 


VI 





1 


0/ 



/o 





1 


1\ 


1 


1 














1 


1 


V 


1 





0/ 


(0 


1 





1\ 


1 





1 





1 





1 





vo 


1 





1/ 



En 



\ I the remaining can be found acting on these five forms with the group of 
symmetries of the square Dg. Combining any four forms, as we have done 
in the previous example, the resulting matrix is always an additive magic 
square, not necessarily normal, and it corresponds to a multiplicative one. 

Theorem 4. fAi, A2, A^, A4 are fundamental forms, or matrices obtained 
from fundamental forms through the action of Dg, then 

8A1 + 4A2 + 2A3 + A4 

is always an additive magic square which corresponds to a multiplicative 
magic square. 
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Given a magic square 8^1 + 4^2 + 2^3 + A4 constructed using our forms, 
we can consider all the permutations of the forms Ai, A2, A3, A4. After a 
permutation, we have still a magic square. Thus we can use the group ^4 
over our forms. 

Example 2. We consider 



/O 1 1\ 
110 
11 

yi 1 0/ 



+ 4 



/I 





1 








1 


1 








1 





o\ 


1 





1 





+ 2 





1 


1 





+ 








1 


1 





1 





1 


1 








1 








1 


1 


Vo 


1 





V 












V 




V 


1 





0/ 



which corresponds to 



/5 

12 
2 



3 
10 
4 
13 



14 
7 
9 




8\ 
1 

15 
6/ 



Permuting the positions of the forms, we obtain 



/I 1 0\ 
11 
11 

yi 1 0/ 

which is equal to 



+ 4 



/lO 
6 
1 

V13 



9 
5 
2 

14 



7 
11 
12 




4\ 
8 
15 
3/ 








1 


A 









1 


o\ 






1 


1 




1 


1 








+ 2 


1 





1 





+ 





1 


1 











1 


1 





1 





1 


1 








1 


V 


1 





0^ 






1 





V 




V 








1/ 



We have seen that using our forms we obtain always an additive magic 
square but it is not necessarily normal. Finally, let us see how to utilize the 
forms in order to generate all and only the normal additive magic squares 
corresponding to normal multiplicative magic squares. The orbit of the 
fundamental forms with respect to the action of Dg are 



A = {Ao,Ai} 

B = {Bq, Bi, B2, B3} 

C = {Co, Ci, C2, C3} 

D = {Do,Di,D2,D3} 

E = {Eq,Ei}. 
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We easily observe that 



/I 



1 
1 
1 
1 



1 
1 
1 
1 



1\ 



Ao + Ai = Eo + Ei=l] 



1 
1 



1 

1 



V 



Furthermore, for any form in the orbits B, C, D there is another form, in 
the same orbit, such that their sum is U. 

We call class the set {A, B,C, D), whose elements are all the magic 
square obtained combining the forms belonging to the orbits A, B, C, D 
(e.g., 8^0 + 4Bi + 2C2 + D3 or 8B2 + 4Ci + 2Dq + Aq are elements of 
(A, i?,C, D)). In the next theorem we show all the classes which provide 
normal additive magic squares. 

Remark 1. We obtain normal additive magic squares, which corresponds 
to all the normal multiplicative magic squares, only from the classes 



We have to made clear some details: 

1. when we choose a form in the orbit C, the forms available in the orbit 
D are only two and vice versa, except for the class {C,C, D, D); 

2. when we have a class with two forms from the same orbits, their sum 
must not be U; 

3. we can not take two times the same form in the same class. 

Considering this remarks, we can count the magic squares obtainable from 
these classes and we check that they are exactly 4^24- 

For the class [A, C, D, E) we can choose 2 form,s from the orbit A, 4 from 
the orbit C, only 2 from the orbit D and 2 from E. Thus from this class we 
can obtain 32 normal additive magic squares, unless of permutations. Thus 



{A,C,D,E) 

{B,B,C,C) 
{B,B,C,D) 
{B,B,D,D) 
{B,C,C,D) 
{B,C,D,D) 
{C,C,D,D) 
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we have 32 ■ 24 = 768 normal additive magic squares from {A, C, D, E), and 
we write 

\{A,C,D,E)\ = 768. 

Similarly we find 



\{B,B,C,C)\ 


= (4 


2 


4 


2) 


6 


= 384 


\{B,B,C,D)\ 


= (4 


2 


4 


2) 


12 


= 768 


\{B,B,D,D)\ 


= (4 


2 


4 


2) 


6 


= 384 


\iB,C,C,D)\ 


= (4 


4 


2 


2) 


12 


= 768 


\{B,C,D,D)\ 


= (4 


2 


4 


2) 


12 


= 768 


\{C,C,D,D)\ 


= (4 


2 


4 


2) 


6 


= 384 



and 



768 + 384 + 768 + 384 + 768 + 768 + 384 = 4224. 



Remark 2. All the magic squares that can be represented through our no- 
tation can be classified and identified by the membership class. 

Remark 3. Such representation allows to construct in a simple way all the 
normal multiplicative magic squares of order 4- 

We conclude this paper with a further example. 
Example 3. Let us consider the factorization 2 • 3 • 5 ■ 67 of 2010. We take 
a magic square in M2010 4 



/2010 5 67 6 \ 

3 134 10 1005 

2 201 15 670 

V335 30 402 1 / 



/2 • 3 • 5 • 67 

3 

2 

V 5-67 



67 



2-3 \ 



2- 67 2-5 3-5-67 

3- 67 3-5 2-5-67 
2-3-5 2-3-67 1 



its image through / is 

/nil 0010 0001 1100\ 

0100 1001 1010 0111 

1000 0101 0110 1011 

\0011 1110 1101 0000/ 



/15 2 1 12\ 

4 9 10 7 

8 5 6 11 

V 3 14 13 / 



This square belongs to the class (C, C,D,D), in fact it can be decomposed 
as follows 

/I 1\ 
110 
10 1 
\0 1 1 0/ 



+ 4 



A 








A 




n 


1 














1 


o\ 


1 








1 


+ 2 








1 


1 


+ 





1 





1 





1 


1 











1 


1 





1 





1 


VO 


1 


1 


V 




V 


1 





V 




^1 





1 


0/ 
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